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We consider a grapliene siieet in the vicinity of a substrate, wiiich contains charged impurities. 
A general analytic theory to describe the statistical properties of voltage fluctuations due to the 
long-range disorder is developed. In particular, we derive a general expression for the probability 
distribution function of voltage fluctuations, which is shown to be non-Gaussian. The voltage 
fluctuations lead to the appearance of randomly distributed density inhomogeneities in the graphene 
plane. We argue that these disorder-induced density fluctuations produce a finite conductivity even 
at a zero gate voltage in accordance with recent experimental observations. We determine the width 
of the minimal conductivity plateau and the typical size of the electron and hole puddles. We also 
propose a simple self-consistent approach to estimate the residual density and the non-universal 
minimal conductivity in the low-density regime. The existence of inhomogeneous random puddles 
of electrons and holes should be a generic feature of all graphene layers at low gate voltages due to 
the invariable presence of charged impurities in the substrate. 

PACS numbers: 81.05.Uw, 72.10.-d, 73.40.-c 



I. INTRODUCTION 

Charge inhomogeneities are known to play an impor- 
tant role in understanding transport properties of semi- 
conductors. The randomly positioned impurity ions give 
rise to a random electrostatic potential and to local mod- 
ulations of the electron density. The corresponding phe- 
nomena have been extensively studied in bulk three- 
dimensional semiconductors and two-dimensional semi- 
conductor heterostructures.'i Recently, it became clear 
that in graphene,^ which is a zero-gap semiconductor, 
the effects of Coulomb disorder on transport are even 
more drastic. 

A particularly interesting property of graphene trans- 
port, dubbed the "minimal conductivity," is a saturation 
of the conductivity, which happens at a low gate voltage 
with a width of AV^ ~ 0.5 eV, which is strongly sam- 
ple dependent. The minimal conductivity value itself is 
controversial with both a universal conductivity^ and a 
sample-dependent non-universal minimal conductivity^ 
being claimed in the experimental literature. Perhaps 
the most remarkable feature of this low-density residual 
conductivity is the approximate saturation of the mini- 
mum conductivity with a voltage width of A , without 
it dropping to zero at = 0, as the density is low- 
ered from an approximately linear voltage dependence of 
the conductivity at high density^i^. Although there are 
transport theories of Dirac fermions in a "white-noise" 
disordered potential- or in clean systems^, which give a 
universal conductivity at the Dirac point, both the mea- 
sured value of the (non-universal) minimal conductivity 
and the observed existence of a saturation region around 
Vg = are in direct conflict with the intrinsic Dirac point 
physics. 

Recent experiments^^ have convincingly proven that 
the observed phenomenon of "minimal conductivity" oc- 



curs due to charge disorder trapped in the Si02 sub- 
strate. Each charged ion produces a potential in the 
graphene layer, which locally mimics a gate voltage and 
thus leads to a non-zero density in the corresponding re- 
gion even if the applied gate voltage is zero (i.e., at the 
"Dirac point" ). If the substrate is electrically neutral, the 
charges in the substrate do produce local density inho- 
mogeneities of electrons and holes, but the total density 
is zero. The physical picture is that of potential valleys 
that define electron puddles and mountains that define 
conducting hole puddles^. For the case of charge neu- 
trality, the point of zero external gate voltage defines the 
percolation threshold, in the sense that the areas where 
the hole and electron densities are non-zero are exactly 
equal to one half of the total area of the sample (which 
is the percolation threshold in two dimensions by dual- 
ity). By changing the external gate voltage, one creates 
an excess of electrons or holes and only one type of car- 
rier percolates. As long as the external gate voltage is 
much smaller than the typical voltage fluctuation due 

to charged impurities, Vg <^ \/v^, the conductivity is 
mostly determined by the random voltage fluctuations 
and is almost independent of Vg. In this case, the con- 
ductivity dependence is expected to have a plateau of 

width \fv^ and be symmetric with respect to — > — Vg. 
The second scenario that should be considered is a sub- 
strate, which is not electro-neutral (which is perhaps a 
more experimentally relevant picture). In this case, the 
dependence of the conductivity on the gate voltage is 
not symmetric with respect to Vg — > — V^, but is shifted 
left or right, depending on the substrate's total charge. 
The percolation picture still holds, and in two dimen- 
sions, the duality between the percolating cluster and 
its shoreline guarantees that the critical gate voltage at 
which the electrons stop percolating is the same as that at 
which the holes begin to percolate. The transition hap- 
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pens at some point Vg ^ V. The width of the minimal 

conductivity plateau is expected to be of order 
At gate voltages much larger than this value, transport 
is expected to be insensitive to the voltage fluctuations 
and is described in terms of the standard transport the- 
ory developed earlier by two of the authors^. We should 
mention that the picture outlined above is entirely classi- 
cal. At the lowest temperatures, the localization physics 
should become important. However, we believe that the 
current experiments are being done in the temperature 
regime where the quantum interference corrections are 
irrelevant - a quantitative agreement between the avail- 
able experimental data at high carrier density and the 
Boltzmann theory of Rcf. 5] strongly support this claim. 

To describe the inhomogeneous density profile and 
transport near the Dirac point in graphene, it is im- 
portant to have a complete description of the statistical 
properties of the corresponding random electrostatic po- 
tential, which creates the inhomogeneous structure in the 
first place. In this paper, we develop a systematic ana- 
lytic theory to derive the relevant probability distribution 
functions of voltage fluctuations and correlators of volt- 
ages due to the charge inhomogeneities. We note here 
that similar potential fluctuations due to randomly dis- 
tributed charges have been considered previously in the 
literature in the context of electron transport in semi- 
conductor heterostructures. E.g., Refs. [lj|9] have stud- 
ied the statistics of potential and density fluctuations 
numerically in the framework of a non-linear screening 
model. Below, we study the voltage fluctuations analyt- 
ically. Using the developed method, we show that the 
PDF of voltage fluctuations is generally non-Gaussian 
and derive explicit analytic expressions for the PDF in 
various limits. Using the obtained analytical results, we 
estimate the typical size of charge-disorder-induced elec- 
tron/hole puddles for a typical graphene sheet. The cor- 
responding results appear to be in a very good agreement 
with the available experimental dataJ^ We also suggest 
a self-consistent (mean-field-like) procedure to estimate 
the typical density at the Dirac point and the remanent 
conductivity near the percolation threshold. We also dis- 
cuss the relation between the suggested percolation-like 
picture of graphene transport near the Dirac point and 
the usual diagrammatic transport theory, which works 
well at high gate voltages. 



II. GENERAL METHOD OF CALCULATING 
STATISTICS OF VOLTAGE FLUCTUATIONS 

Let us consider a two-dimensional plane in the vicinity 
of a substrate containing randomly positioned charged 
impurities. Carriers in the 2D plane screen the bare 
Coulomb potential produced by the charges in the sub- 
strate. The screening properties and the corresponding 
effective potential depend on the nature and the density 
of the carriers, but for the purpose of this section, the ex- 
act form of the potential is not important. Let us denote 



it as v{r). 

To derive the PDF of voltage fluctuations, we exploit 
a method, used e.g. by Larkin et alJ^ in the context of a 
mean- field approach to random spin systems. First, we 
define the PDF as follows 



P[V] 



(1) 



where V is the random voltage, the index i in the sum 
labels the charges, v{pi) is the potential produced by the 
i-th impurity, and the angular brackets correspond to 
averaging over all possible positions of impurities. Using 
the standard representation of the (^-function, we get 



P[V] = J ^e^""' {J[e^pHXv{p)] 



(2) 



At this point, we assume that the impurities are uncorre- 
lated, which allows us to exchange the averaging opera- 
tion with the product over the impurities in Eq. We 
obtain 



P[V] 
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exp [-iXv{p)]) 
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where A^imp is the number of impurities in the substrate. 
Now, we average over all possible positions of impurities, 
which have the equal probabilities. This implies sim- 
ply evaluating an integral over the entire two-dimensional 
area A. 

{exp[-iXv{p)]) = 1 - ly -exp [-iXv{p)]} , 

A 

(4) 

where we performed a trivial operation of subtracting 
and adding a unity. Using Eq. (jl]), we see that in ther- 
modynamic limit, the PDF can be written as 



P[V] 
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/ d p{l - exp [iXv{p)]} 



,(5) 



where nimp = Ni^p/A is the concentration of impuri- 
ties. Note that this result does not depend on the spe- 
cific form of the potential and is completely general. We 
emphasize here that Eq. ([5|) describes a random poten- 
tial due to a charged substrate, which contains disor- 
der of just one type. If the substrate contains charges 
of different types, at the level of Eq. (|4]), we have to 
perform another averaging over the distribution function 
of charges, C[Q]. In particular, if the substrate is elec- 
troncutral and contains charges of two types of opposite 
signs itQo, then the charge distribution function is sim- 
ply C[Q] = \ {8{Q - Qo) + 6iQ + Qo)], which leads to 
the following PDF for the electrically neutral substrate: 



Po[V] - 
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- cos [Xv{p)]} 



.(6) 
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To describe the distribution and sizes of disorder- 
induced droplets, we also consider another PDF function, 
which determines the distribution of voltages in two dif- 
ferent points in the film. Similarly to Eq. ([1]), we define 
it as 



J 



(7) 

where pij = pi — Pj is the distance between the two 
points. This correlation function ([7]) determines the con- 
ditional probability of finding a voltage V2 in the point 
P12, provided that in the origin, the voltage is equal to 
Vi. Thus, the decay of this correlation function will de- 
termine the typical size of the charge disorder-induced 
puddles of electrons or holes. 

Following the simple steps, which led to Eq. ([5]), we 
obtain for Eq. ^ 



on^ ^r 1 f d\id\2 



(8) 



"imp / d?Pa {1 - exp [i\2v{p2o) - iXiv{pai)\} 



X exp 



Formally, Eqs. ([5]) and ^ completely describe the 
probability distribution of voltages for a given effective 
interaction, v{v). However, exact analytical calculation 
of the full PDFs is possible only in the simplest cases of 
pure Coulomb interaction and large-distance asymptote 
of the linearly screened Coulomb (see below Sec. IIIII) . 
To correctly interpolate between the two behaviors, one 
has to consider the general form of the screened inter- 
action potential. The latter has a relatively complicated 
structure in real space and generally it is not possible to 
calculate the full PDF analytically. However, it is possi- 
ble to evaluate exactly the relevant moments of the PDF. 
These moments can be calculated from Eqs. ([5]) and ([5]), 
by simply evaluating derivatives with respect to the aux- 
iliary parameter A, 
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In a similar way, one can determine the real space corre- 
lation function of the random voltages 

[V{v)-V\ [V{Q)-V\ =ni^p / ^i\q)e^^-\ (10) 
which is valid for any effective potential v{q). 



III. DERIVATION OF PDF FOR SPECIFIC 
MODEL INTERACTIONS 

In this section, we demonstrate the application of the 
method of Sec. |TT] by considering three specific model 
forms of the interaction v{p). 



A. PDF of voltage fluctuations for bare Coulomb 
potential 

We start with the bare Coulomb interaction vq{p) = 
Qo/ p, where ±Qo are the charges of the impurities. We 
have to evaluate the following integral [see Eq. dS])]: 



/o(A) = 2m._ 
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Apparently, this integral diverges at large distances. 
So, we expand the cosine term and get /o(A) = 

Pmax 

TrrtimpA^Q^ / dp/p = TrnimpA^Qg In (pmax/Pmin)- To 

Pmin 

regularize the divergences in the integral, we note that 
the minimal distance can not be smaller than the distance 
from the substrate. On the other hand, the large-distance 
divergence is connected with the long-range nature of the 
Coulomb interaction and is regularized by screening char- 
acterized by a Thomas-Fermi screening wave- vector qg. 
Summarizing, we write the PDF corresponding to the 
pure Coulomb interaction 



Pbarc[^] - / ^ exp \iVX - a(A)A2] , 



(11) 



with 



a(A) 



^imp 



Q^ln [q^ ^max {d, |A|Qo}] 



Note that this coefficient depends on A very weakly. We 

can substitute A with a typical A ~ max |y , 

After this simplification the integral (fTTj) can be calcu- 
lated and we find a Gaussian PDF of voltage fluctuations 
with the variance V"^ = 2a. We see that the typical volt- 
age is enhanced as compared with the nai've estimate, 

1 /2 

"-impQoi by a large logarithmic factor. 



B. PDF of voltage fluctuations for large-distance 
asymptote of screened interaction in 2D 

Another model interaction, which allows analytic 
treatment is the large-distance asymptote of the screened 
interaction. Due to a wek screening in two dimensions, 
it decays as a power law 

VscAp) = 'iQo/q^^p^. 

To determine the statistics of the random voltage, we 
need to evaluate the following integral [see Eq. ^] 



^cr(A) = 27rn 



00 

imp J dpp 





cos ^— - 



This integral does not have any infrared problems due 
to screening. After some algebra, we find the following 
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PscriV] = ^/ dxcoa (^^x^ e-^"'\ (12) 



where (3 = ni,„p[r(l/3)7r/2] [AQo/qiy'^^ is a parameter, 
which depends on the carrier density and r(z) is the 
Gamma function. From Eq. (fT^ . one can see that in the 
case of a screened interaction the PDF of voltage fluctu- 
ations is strongly non- Gaussian. One can check that aU 
even moments of this PDF diverge (all odd moments are 
trivially zero due to assumed electroneutrality) , which 
corresponds to a singular behavior of the model poten- 
tial at small distances. Obviously, this divergence is an 
artifact of our approximation about the interaction po- 
tential, which is justified only at large distances. 



IV. SCREENED COULOMB POTENTIAL 
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FIG. 1: (Color online) This figure shows the representative 
probability distribution functions of random voltages aris- 
ing from the screened potential (|13|l of randomly distributed 
charges in an electrically neutral substrate for three different 
densities (other parameters are fixed). 



To correctly interpolate between the small distance 
(large-y) asymptotics of Sec. IIII Al and large-distance 
(small-y) behavior of Sec. IIII B( we need to consider a 
general form of the effective interaction, which includes 
both screening and a finite width of the spacer layer. Be- 
low, we specialize to the case of linear screening, where 
we use the following form of the effective potential 



v{q) 



2nQie 



— qd 



9 + 9s 



(13) 



where is the Thomas- Fermi screening length and d 
is the distance between the 2D graphene plane and the 
substrate. Note that the screening length is graphene 
is density-dependent. If there is no external gate volt- 
age, there is no density to start with. However, it in- 
evitably appears due to the random voltage fluctuations 
and can be determined self-consistently in the framework 
of a mean-field-like theory (see Sec. Ulland Ref. [IH). In 
this section, we simply take the quantity qs as a free pa- 
rameter of the theory. 

Due to a complicated structure of the potential in real 
space, it is not feasible to calculate the full PDF (O ana- 
lytically. However, it is easy to perform a numerical eval- 
uation of the corresponding integrals. A result of such 
calculation is shown in Fig.[TJ where representative PDFs 
are plotted for three different densities. We note that 
different portions of the curves can be very well-fitted 
with Gaussians, but the entire dependence is definitely 
more complicated than a simple Gaussian. Moreover, if 
y ^ (a charged substrate), the distribution function is 
not even symmetric with respect to 5V = {V — V^, i.e. 
P[5V] ^ P[-5V]. 

Even though the full PDF can not be determined an- 
alytically, one can evaluate explicitly its moments using 
Eqs. ^ and PH)) . In particular, we calculate the sec- 
ond moment for the case of Thomas-Fermi screening, 



6V^ = 



= 27rni,„i 



the result is 



Ql [El (2g,d)e2«='' (1 + 2qJ) - l] , (14) 

where Ei{z) — t^^e^*^dt is the exponential integral 
function. Obviously, the statistical properties of the po- 
tential strongly depend on the parameter qsd. Below we 
give asymptotic expressions for SV^ in the limits of large 
and small qsd 



(2g,d)" 



In 



if qsd > 1; 
if qsd < 1, 



(15) 



where 7 « 1.781 is the Euler's constant, e w 2.718 is 
the base of the natural logarithm. We note that the 
last expression essentially reproduces the result obtained 
above for the bare Coulomb potential [see Eq. ((TT|) and 
the text below it]. 

Next we study the correlation function {V{r)V{Q)), 
which is related to the typical size of electron/hole pud- 
dles. To simplify the notation we define a new function 
C(r) according to the following relation {V{r)V{0)) — 

V + 27rninip(5§C(r). Using Eq. pO|) . one can calculate 
analytically all coefficients in the Taylor series for the 
correlator C(r) 



C(r) = Co(2g,d)+^(-l)™(9,r) 

m>0 



\2m 



{2m - 1) 



-an(2qsd), 
(16) 

where C™(x) = [I / {2my]{dl"'+^)[xe'= Ei{x)]. Note that 
Co reproduces the expression in the square brackets in 
Eq. (fH|) . In the limit qsr ^ 1, the correlator behaves 
as C{r) « Co — (Ci/2)(gsr)^. In the opposite limit of 
qsT ^ 1, one finds that the correlation function decays 
as a power law C(r) 2|1 — qsd\{qsr)~^ . The general 
form of the correlator has been calculated numerically in 
Sec. IV Al and is shown in FigO 
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APPLICATION TO GRAPHENE 2 



A. Size of the puddles 



Now, we apply the general results obtained above 
specifically to the case of a graphene sheet in the ex- 
perimentally relevant parameter regime. We estimate 
the density of charged impurities, the parameter qgd, 
which determines the behavior of the correlation func- 
tion C(r), the typical size of a fluctuation- induced elec- 
tron/hole droplet due to the charge impurity fluctua- 
tions, and the typical width of the region where the con- 
ductivity is expected to be almost a constant. We esti- 
mate these quantities for typical graphene samples with 
mobilities /i ~ W^crn^ /Vs. The deviation from the semi- 
classical Drude-Boltzmann theory occurs at carrier den- 
sities n ~ nimp « 5 X lO^^cm"^ This leads to the 
density of states vo ~ 10"'^''cm"^meV~^ and the Fermi 
energy 80 meV. Assuming that the distance be- 

tween the graphene plane and the substrate is c? ^ Inm, 
we find Qgd ~ 0.5. The corresponding correlation func- 
tion is plotted in Fig. 2. This correlation function de- 
termines the size of a disorder-induced puddle, which 
can be roughly estimated as I ~ lOnm. We also can 
determine the average value of the random potential as 

V = riimp/i^o ~ 0.05 eV and VW^ « 0.1 eV. The for- 
mer number determines the asymmetry in the conduc- 
tivity dependence, cr(V^), while the latter estimates the 
width AVg of the constant conductivity plateau. These 
numbers are in a reasonable agreement with the available 
experimental data and validate the basic model that the 
minimum conductivity in graphene is connected with the 
density inhomogeneities caused by the charged impurities 
in the Si02 substrate. We note two direct experimen- 
tal consequences of our theoretical findings: (1) Dirtier 
samples should have larger vales of AVg] implying higher 
(lower) mobility samples should have smaller (larger) val- 
ues of AV^. (2) In addition, AVg has a nontrivial depen- 
dence on d, the impurity location [see, Eq. (|15p ]. which 
can, in principle, be experimentally verified. 



While the voltage (or density) width of the residual 
minimal conductivity regime is explained naturally in our 
theory of impurity induced inhomogeneous electron-hole 
puddle picture, getting the exact values of the residual 
density and the minimal conductivity will require a de- 
tailed calculation involving percolation through a net- 
work of random 2D electron and hole puddles i^^i^^ The 
following section provides estimates of these quantities in 
the framework of a self-consistent "mean-field" approach, 
which should be quantitatively reliable away from the 
precise percolation point. 
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FIG. 2; (Color online) The correlation function C(r) [see 
Eq. (|10|) and text for the definition] for the charge impurity 
model with linear Thomas-Fermi screening. The solid line 
shows numerical integration for Qsd = 0.5, while the dashed 
lines correspond to the asymptotic analytical expressions for 
small and large gsr. The chosen value of qsd = 0.5 is deter- 
mined by the density of carriers, where the experimentally 
observed behavior deviates from the predictions of the classi- 
cal Drude-Boltzmann theory. 



B. Mean-field approach and relation to 
diagrammatics 

In this section, we briefly outline possible approaches 
to relate the voltage fluctuations to observables and, 
in particular, discuss the regime of high gate voltages, 
where the usual diagrammatic perturbation theory be- 
comes quantitatively applicable and show a possible re- 
lation between the discussed charge disorder fluctuation 
phenomena and the diagrammatic approach. We also 
formulate here a self-consistent mean- field- like approxi- 
mation, which allows us to roughly estimate some observ- 
able parameters in the regime where perturbation theory 
breaks down. 

While it is intuitively clear that the local voltage 
fiuctuations lead to a non-vanishing density in the 
two-dimensional plane, the exact relation between the 
voltage distribution and the density is non-trivial. 
Generally, one has to solve the many-body Schrodinger 
equation taking into account the long-range disorder 
potential and Coulomb interaction between carriers, 
which is an extremely complicated if not impossible task 
to accomplish. Efros et alii have previously formulated 
a phenomenological approach in which the density is 
determined from minimization of the energy of the sys- 
tem written in terms of the position-dependent density 
n = J d'^r{^V{r) [n(r) -n\+ E[n{r)]}, where the first 
term describes the Coulomb interaction between the 
impurity charges and the induced carriers and the local 
density, while the second term represents the energy 
of a uniform liquid (which includes the kinetic and 
exchange correlation effects). Formally, this approach 
would allow to interpolate between the regimes of large 
and small carrier density; but in the former case, i.e., 
if the chemical potential is large, the effect of charge 
disorder fluctuations is a small perturbation and the 
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FIG. 3: Fig. 3a is the pictorial representation of the 
disorder-induced density fluctuations {5n{v)5n{v')) , where 
the dashed line represents the Coulomb impurity potential 
and the shaded vertex is a ladder of these impurity lines. 
Fig. 3b is a diagram, which contributes to the correlator 
((5n(ri, r2)5n(r3, r4)) of the polarization operators, which de- 
termine the screening properties and the effective interaction 
(see text). 



aforementioned phenomenological approach should 
reduce to the usual microscopic perturbation theory. 
E.g., one can determine the relevant density-density cor- 
relation function {5n{r)6n{v'))) by considering diagrams 
in Fig. 3a, where the dashed line represents the electron 
scattering off of the Coulomb impurity potential and the 
shaded vertices imply a ladder of such charge impurity 
lines. We note that to avoid infrared divergencies, 
one has to consider screening of the impurity potential 
by carriers, C/(r,r') ~ v{y — r') + * 11 * C/](r, r') 
(where v is the bare potential and 11 is a polarization 
operator). Strictly speaking, the screening properties 
and the interaction between two charges are random 
quantities themselves, which depend on the positions 
of the charges. In particular, the polarization operator, 
which determines the screening properties is a random 
matrix, n(r, r'), which depends on two coordinates. 
The probability distribution function of this matrix (in 
the Gaussian approximation) is a functional i^l^II] = 

exp [-(1/2) 2,3,4 '^n(l, 2)C-i(l, 2; 3, 4)<5n(l, 2) , 

where C~^(l,2;3,4) is the inverse correlator 
((5n(ri, r2)(5n(r3, r4)), which is described by the di- 
agrams in Fig. 3b. Clearly in the high-density regime, 
the average screening length q'^^ will be determined 
just by the average density, with small "mesoscopic" 
corrections. Similarly in this regime, where the Fermi 
wave-length is the smallest length scale in the prob- 
lem, the local density in a particular point will be 
determined just by the local voltage which can be 
viewed as a local chemical potential in this point. 
Therefore, the average density fluctuation (up to a 
numerical constant of order one) should be of the order 

[where ?i(/i) is the 
dependence of the density on the chemical potential in 
the uniform case and SV'^ is given by (fT4|) ]. 

When the external gate voltage and thus the chemi- 
cal potential and the Fermi-momentum decrease, the mi- 
croscopic perturbation theory of density fluctuations in 



FIG. 4: The first few diagrams in the perturbative expan- 
sion of the conductivity. The dashed lines correspond to the 
scattering off of Coulomb disorder potential. 



graphene breaks down. The Green's functions and the 
screening properties of the impurity potential in Fig. 3 
both depend on the average density, which vanishes at the 
Dirac point and the corresponding computation becomes 
ill-defined. However, it is possible to formulate a reason- 
able and simple self-consistent scheme, which would al- 
low to get an order-of-magnitude estimate of the density 
in the vicinity of the Dirac point. Namely, one assumes 
that the Green's functions and the screening length cor- 
respond to a uniform density ^ \/ S^n. Then, using these 
objects one can calculate the density fluctuation, which, 
as argued above, should be of order 

\/ S'^n n{ii ' 

But in this equation, the right hand-side of this relation 
itself depends on S'^n via the screening length [see, 
Eq. (fT4|) of Sec. lIV| . This leads to an algebraic non-linear 
self-consistency equation for the typical density fluctua- 
tion. The corresponding calculation has been done by us 
and is described in detail in Ref. [11] (the corresponding 
estimate for the residual density is n ~ 0.1 — Inimp). 

A more important experimentally measurable quan- 
tity is the conductivity. At very high gate voltages, the 
usual semiclassical Drude-Boltzmann theory works per- 
fectly well. As the chemical potential is lowered, the 
"mesoscopic fluctuation" effects start to appear. In the 
regime, when these effects are a small correction to the 
classical conductivity, one should be able to obtain them 
in the framework of the conventional diagrammatic tech- 
nique. This implies going beyond the Drude approxima- 
tion (no crossed impurity lines) and calculating "quan- 
tum corrections" to the conductivity. A formal calcula- 
tion would select the maximally crossed weak-localization 
diagram, which diverges in the zero-temperature limit. 
However, the flattening of the linear-in-T^ conductivity 
most likely has nothing to do with the weak localization 
or any other quantum interference effect of this type. The 
effect of minimum conductivity is "quantum" only in the 
sense that diagrammatically it appears beyond the clas- 
sical Drude approximation and is due to diagrams with 
some crossed impurity lines. We note that in the case 
of long-range Coulomb disorder, there are other dimen- 
sionless parameters apart from kpl (namely, which for- 
mally should be assumed small for perturbation theory to 
hold). If the dephasing time is short enough, other types 
of diagrams (e.g., impurity vertex corrections) may be 
more important than the weak localization diagram. We 
believe that in the recent experiments, which observed 
the minimum conductivity, the latter scenario is real- 
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ized. The fate of the weak locahzation correction is a 
very interesting but separate question, which seems to be 
irrelevant to the effect of the flattening of conductivity 
at low gate voltages. The initial deviation of the (j{Vg)- 
dependence from the straight line can be obtained from 
the usual perturbation theory by considering just the 
next-to-leading diagrams with a finite number of crossed 
Coulomb impurity lines (see Figs. 4b and 4c). The cor- 
responding calculation will be published elsewhere. 

As noted above, in the limit of zero or very small gate 
voltage (in the very vicinity of the Dirac point), the 
microscopic perturbation theory breaks down and only 
qualitative descriptions are available. In this regime, 
transport occurs via a percolation network of electron 
and hole puddles. Cheianov et alJ^ recently proposed a 
random resistor network model and calculated a scaling 
behavior of the conductivity near the percolation thresh- 
old. While the corresponding scaling function may be 
universal, the value of the minimal conductivity itself is 
non-universal (e.g., it depends on the choice of the un- 
derlying lattice) and can not be determined within a ran- 
dom resistance network model. The possibility that dis- 
order associated with random charged impurities could 
lead to a non-universal graphene minimum conductivity 
was pointed out in Ref. [5| , and the value of the minimal 
conductivity was estimated in Ref. flTl , by assuming that 
it is of the order of Drude-Boltzmann conductivity with 
the density determined by the value of the typical den- 
sity obtained via the self-consistent procedure explained 
above. We note that this estimate agrees very well with 
the available experimental data-^ (e.g., both the high- 
voltage Drude behavior and the minimal conductivity 
can be determined using just one fitting parameter - the 
density of impurities in the substrate). 



VI. CONCLUSION 

In this paper, we developed a simple and general 
method of calculating statistical properties of voltage 
fluctuations due to charge disorder in the vicinity of a 
conducting two-dimensional plane. We used these results 
to describe the low-density properties of a graphene sheet 
near a substrate with Coulomb impurities. In particular, 
we estimated the typical size of the electron and hole pud- 
dles and the width of the minimum conductivity plateau. 
We also proposed a self-consistent approximation scheme 
to estimate the typical residual density and conductivity 
at the Dirac point. The corresponding results appear to 
be in a good agreement with experimental data.— 

While our theory based on density fluctuations by ex- 
trinsic charged impurities invariably present in the sub- 
strate provides a reasonable mean-field description of 
the observed graphene non-universal conductivity min- 
imum plateau around the charge neutrality point, an im- 
portant open question remains on the transport behav- 
ior precisely at the percolation critical point Vg = Vc 
where infinite electron and hole percolation paths open 



for transport throughout the system. Neither the mean- 
field theory of Ref. [ll| nor the perturbative diagram- 
matic approach can qualitatively access the critical point 
where the density fluctuations are very large. Unfortu- 
nately, the problem of determining the minimum con- 
ductivity at Vc remains inaccessible to the theory. We 
should mention that recently Cheianov et al.^^ studied 
the critical percolation problem by effectively including 
tunneling between the neighboring electron/hole puddles 
within a lattice random-resistance model, however their 
numerical analysis has not led to any estimate of the non- 
universal non-universal conductivity at the percolation 
critical point. Experiment seems to indicate a smooth 
and continuous behavior as one goes from the Vg < Vc to 
the Vg > Vc regime with the minimum conductivity re- 
maining approximately a constant over a finite interval of 
voltage 5Vg, i.e. in the regime Vg — 5Vg <Vc <Vg + 5Vg, 
as the voltage swept through the charge neutrality point. 
The "mean-field theory" seems to provide a reasonable 
quantitative estimate of both (7,nin ^ '^(^c) and the value 
of 5Vg itself as well as the conductivity cr (|V^|) \Vc\ 
far from the electron-hole puddle percolation regime. 
Whether there is some remnant universal quantum con- 
ductivity behavior, arising from the special quantum- 
mechanical properties of the chiral Dirac-like graphene 
band dispersion, remains an important open question. 
The smooth and continuous experimental behavior of 
graphene transport through the charge neutrality point 
empirically argues against any universal behavior, but 
we can not rule it out on the current theoretical grounds. 
In fact, closely connected with the critical percolation 
transport is the question of (the absence of any observed) 
Anderson localization in graphene. The strong disorder 
associated with the random charged impurities should 
lead to strong Anderson localization in graphene. Ex- 
perimentally, however, no such strong localization effect 
has ever been seen in graphene and even the most disor- 
dered graphene samples seem to have a finite minimum 
conductivity which is relatively temperature independent 
around T < 100 mK. Why this is so is not understood 
theoretically. It is possible that the crossover to strong 
localization occurs at very low temperatures but more ex- 
perimental and theoretical work is needed to settle this 
important question. The issues of the critical percola- 
tion transport near the charge neutrality point and the 
apparent absence of Anderson localization in graphene 
are both beyond the scope of this work. 

An important conclusion, which follows from our work 
is that there will be no percolation induced metal- 
insulator transition in graphene as there is in 2D semi- 
conductor systemsi^ and the graphene layer will always 
conduct at low density (leading to a residual conductiv- 
ity), as long as there is no Anderson localization, simply 
because the electron and hole percolation densities are 
exactly equal by duality in two dimensions, and there- 
fore, either the electron or the hole channel always con- 
ducts independent of how strong the impurity-induced 
inhomogeneous voltage fluctuations are. This residual 
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conductivity is not universal and depends on the density 
of charge impurities. 

FinaUy, we emphasize that the theoretical technique 
we develop in this paper for calculating the non- 
Gaussian probability distribution function of charged im- 
purity induced voltage fluctuations is completely gen- 
eral, and can be used for other 2D systems, e.g. 
modulation-doped semiconductor heterostructures and 
quantum wells, where such fluctuation-induced density 
inhomogeneity and puddle formation effects^^ are known 
to be important, often leading to percolation- induced 



conductor-to-insulator transition. We believe that the 
theoretical framework developed in this work can be 
adapted to 2D disordered semiconductors providing in- 
sight into the nature of transport in these systems. 

Note added: After the first version of this manuscript 
had been submitted, we received a preprintiS, which 
reported an explicit experimental observation of the 
electron-hole puddles in graphene and provided an es- 
timate of a puddle size in quantitative agreement with 
our theory. Acknowledgements: This work is supported 
by US-ONR. 
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